Angular distributions of atomic photoelectrons produced in the UV and XUV regimes 



S. Bauch and M. BonitzQ 

Institut fur Theoretische Physik und Astrophysik 
Christian- Albrechts-Universitat zu Kiel, D-24098 Kiel, Germany 
(Dated: June 27, 2008) 

We present angular distributions of photoelectrons of atomic model systems excited by intense 
linearly polarized laser pulses in the VUV- and XUV-regime. We solve the multi-dimensional time- 
dependent Schrodinger equation for one particle on large spatial grids and investigate the direction 
dependence of the ionized electrons for isotropic s-states as well as p-states. Although the pondero- 
motive potential is small compared to the binding energy of the initially bound electron and the 
photon energy of the exciting laser field, richly structured photoelectron angular distributions are 
found which sensitively depend on the laser frequency and intensity as well as on the number of 
absorbed photons. The occuring shapes are explained in terms of scattering mechanisms. 
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I. INTRODUCTION 

Photoemission experiments are an important tool for 
the investigation of electronic properties of matter from 
single atoms and molecules to condensed matter systems, 
see e.g. [T] for an overview. 

The first experimental investigation of photoelectron 
angular distributions (PADs) from atoms in the gas phase 
using ultra-violet light was published in 1930 |2J. The 
shape of the PAD followed the cos 2 dipole-like shape pre- 
dicted by early quantum mechanics [3] which stated that 
the most favorable direction of emitted electrons is along 
the field polarization axis. 

The technical development of new light sources, in the 
first place the invention of the laser [4], with the abil- 
ity of creating high-intensity and monochromatic pulses 
gave new momentum to this topic, giving access to non- 
linear processes. An (ionization) process is called non- 
linear if the photon energy huj is smaller than the ion- 
ization potential I p of the atomic state considered. Ac- 
cording to the classification by Keldysh [5] two main pro- 
cesses can be distinguished. Depending on the parameter 
7 = Ip/2U P (where U p = e^E^j '4m e uj 2 is the pondero- 
motivc energy of the electron in the changing electric field 
of the laser with amplitude Eq) the electron is freed into 
a continuum state by tunnel ionization (7 < 1) or by 
multi-photon (MP) ionization (7 > 1). 

After its ionization the electron can absorb additional 
photons if the intensity of the laser is sufficiently high. 
The first experimental evidence of such above-threshold 
ionization (ATI) by one additional photon in intense 
infrared (IR) laser pulses has been observed in the late 
70's In this case, the photoelectron is ionized by a 
MP or tunneling process and absorbs additional photons 
from the light field to increase its kinetic energy leading 
to a peak-like structure in the electron spectra, where 
each peak is separated by the photon energy hu of the 
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exciting laser. Soon, better laser systems and higher 
accuracy in experiments showed the expected sequence 
of peaks in the continuum [71 [5] which form well-known 
plateau structures. The quasistatic or two-step model 
[H Uni E] explains the underlying dynamics. First, the 
electron is ejected into a continuum state via tunneling 
or MP processes with nearly no kinetic energy, being 
subsequently accelerated by the electrical field of the 
laser. The main structures of ATI spectra can be 
explained by rescattering processes where parts of the 
electronic wave function are driven back to the parent 
ion. This effect is the origin of such famous and actively 
studied processes as high harmonics generation (HHG), 
non-sequential double (or multi) ionization and the 
formation of the characteristic cut-offs in the energy 
distribution of ATI electrons. 

A large amount of information regarding electron dy- 
namics in atoms and ionization processes is obtained 
by analyzing the corresponding angular distributions of 
electrons. First investigations of multi-photon PADs re- 
vealed intensity-dependent structures |12j . According to 
the simple quasistatic model the general thought was that 
the PADs are more and more peaked along the polariza- 
tion axis of the laser field with increasing photon order 
of the process, which leads to occupation of final states 
with higher angular momentum. But the experimental 
situation soon changed with the appearence of additional 
structures referred to as side lobes\l'S\, jets and wings [T3] 
which appear at characteristic regions in the ATI photo- 
electron spectra, depending on the ponderomotive en- 
ergy. The angular distributions have been widely used to 
identify the involved high- lying Rydberg states ("chan- 
nel switching" ) whose angular momentum significantly 
affects the final angular momentum state of the photo- 
electron, and such the PAD [T5] . 

Much work has been performed in the regime of in- 
tense (< 10 15 W/cm 2 ) IR pulses. A recent comparison 
between experiment and theory showed perfect agree- 
ment |16j . On the other hand, multiphoton ionization 
of rare gas atoms in the regime of high photon energies, 
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which allows for direct ionization of one electron by a 
single photon, has been investigated at the free electron 
laser (FEL) facility at DESY in Hamburg which gave 
the first experimental evidence of a multi-photon process 
in the high photon energy regime |17l 118] , Up to now, 
there still exist discrepancies between theory and exper- 
iment [19 which may also be attributed to still limited 
experimental work in the field since the high-intensity 
VUV and XUV sources are presently under construction. 
Many predictions for further investigations already exist 
[201 EH 123 [231 [2H HSl HSl HZl HE] but, to our knowledge, 
no information on the expected ADs of the photoclcc- 
trons is available. Since the ponderomotive potential is 
of the order of meV in this regime, ponderomotive scat- 
tering should play only a minor role. Additionally, the 
ionization happens directly into the continuum - no chan- 
nel switching due to ponderomotive shifting of resonant 
Rydberg states, affects the final electron state. Hence, 
the ADs should be dominated by the initial state of the 
atom and, with increasing laser intensity (and therefore 
higher electrical field strength), by rescattering effects 
with the parent ion. 

The aim of this paper is the investigation of (multipho- 
ton) PADSs at high photon energies. It is organized as 
follows: In Section [II] we explain our general theoretical 
approach to obtain angular distributions of photoclcc- 
trons which is followed by an explanation of the chosen 
model system in Section |III| Our results are presented 
in Section [IV| in detail for various laser parameters. The 
work ends with a discussion of the effects which con- 
tribute to the PADs and gives a physical explanation of 
the structures in terms of scattering processes by utilizing 
a simple analytical fit formula. 



II. METHOD 
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FIG. 1: Schematic drawing of the used detector implemen- 
tation. The wave function is saved at a large fixed radius r 
after which it is damped by an imaginary absorbing potential. 
Typical detector parameters are r > 25.0 a.u and more than 
1500 detector points. The atom is located in the center region 
at r = 0. 
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FIG. 2: (color online) Density of the initial states of the regu- 
larized Coulomb potential, Eq. |3| with k = 0.1, in the x — y- 
plane calculated via TDSE propagation in imaginary time 
direction. The lower numbers give the energy eigenvalue of 
the state. The nomenclature of the states is chosen in analogy 
to the orbitals of the hydrogen atom. 



We solve the two-dimensional time-dependent 
Schrodinger equation (TDSE) [25] 

.(1) 

for one particle on large spatial grids by means of an im- 
plementation of the Crank- Nicolson procedure in combi- 
nation with the operator-splitting technique which gives 
access to the solution of the two-dimensional problem. 
The potential term V = V atora + V int is given by a time- 
independent part V a tom(^, y) and a time-dependent inter- 
action part Vi n t(x, y, t) which describes the external laser 
field. 

The electro-magnetic field is treated classically within 
the dipo/e-approximation. Then, the interacting part of 
the Hamiltonian for a linearly polarized laser field along 
x-direction reads in length gauge 

Mntfov,*) =exp f- (t ~2 0) ) E xcos(u{t-t )) , (2) 



where the envelope is assumed to be of Gaussian shape, 
r describes the pulse duration (full width at half max- 
imum FWHM), to the time of maximum electrical field 
strength Eq and u> is the photon energy. 
In order to minimize the computational grid sizes a dissi- 
pative absorbing potential with an imaginary part is in- 
cluded at the boundary of the grid. The decrease of nor- 
malization of the wave function due to absorption gives 
then an estimate of the total ionization rate. Special at- 
tention is paid to the calculation of energy spectra of the 
ionized electrons. We implemented a realistic detector- 
like setup (cf. Fig. [lj which allows for the calculation of 
energy- and angle-resolved photoemission spectra. The 
two-dimensional wave function ^(t, r — \J x 1 + y 2 ) is 
saved at the fixed radius r, where r has to be chosen 
large enough to avoid near-field effects, for all time steps 
and hence the energy spectrum on the detector can be 
obtained by a Fourier transform with respect to the time 
t. After its detection the wave function is damped by 
the imaginary part of the potential, cf. Fig. [l] We used 
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FIG. 3: (color online) Angle- and energy-resolved spectrum 
of the ionization of the 2s-state for two photon energies: (a) 
lj — 0.5 a.u. and (b) ui = 1.0 a.u. The electrical field strength 
is fixed at Eo = 0.1 a.u. corresponding to a laser intensity of 
I = 3.5 ■ 10 14 W/cm 2 . The direction of the linearly polarized 
laser field is indicated by the (red) dashed line. 



r > 25.0 and approximately 1500 detector points allow- 
ing for a high angular resolution. During the data pro- 
cessing, the resolution is reduced to an angle element of 
= 2tt/256. 

As initial conditions the cigenstates of the considered 
potential are used. They are obtained via propagation 
of the TDSE in imaginary time direction (time step iAt) 
[30] which is a procedure similar to self-consistent Hartree 
Fock calculations. Higher-lying states (above ground 
state) are constructed by an additional orthogonalization 
procedure at each imaginary time step. 



III. MODEL SYSTEM 

The binding potential of the atom is represented by a 
two-dimensional regularized Coulomb potential 



V atom (x,y) 



Z 



\J x 2 + y 2 + k 2 



(3) 



in our simulations, n is a small cut-off parameter pre- 
venting the singularity at x = y = 0. This is necessary 
for the numerical treatment on the chosen cartesian 
grid. Furtheron, we investigate a hydrogen-like atom 
[Z = —1.0 a.u. ) with a chosen regularization of k = 0.1. 
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FIG. 4: (color online) The same as Fig. [3] but for (a) the 
p x - and (b) the p^-state, excited with a laser frequency of 
lo = 0.5 a.u. 



Throughout this paper the ionization process of the 2s 
and the 2pa;-like and 2p y -like states (cf. Fig. [2l will be 

ncalfi 



investigated. These chosen orbitals are typical for all s 
and p states. The orientation of the states with respect 
to the laser polarization is generated by the imaginary 
time propagation and is - in some sense - arbitrary. 
The explicit alignment of the states will be discussed 



lateron, see Section IV B The deeply-bound Is ground 
state - being present in the construction of eigenstates 
via the imaginary time propagation of the TDSE - 
is non-physical due to the smoothing of the potential 
and is not considered for ionization in the following 
(l-EWdl > 3 a.u.). Also its binding energy exceeds the 
photon energy several times, hence no or only small 
fractions could be ionized (the intensity of most of the 
laser fields considered is too small to allow for significant 
multiphoton ionization from this state). 

The introduced artificial screening by the parameter k 
lifts the degeneracy of the 2p-states and the 2s-state, as 
can be seen in the printed eigenenergies in Fig. [2] In anal- 
ogy to the full hydrogen atom we will call the constructed 
states with / = 1 , for simplicity, p^-state and p y -state, al- 
though they are not aligned along the corresponding axis. 



In Section IV B the dependence of the angular distribu- 
tion of photoelectrons on the explicit alignment of these 
states with respect to the laser polarization axis will be 
investigated in detail. 
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IV. RESULTS 

A. Angular dependence of different photon lines 
for s-states 

We consider pulses of duration r = 30 fs to closely 
reproduce typical experimental situations with different 
photon energies and intensities. Since we use UV photon 
energies lu in our calculations the chosen duration should 
not affect the results significantly as many oscillations 
of the electrical field are present and we are far of the 
regime off few-cycle pulses. 

Figs. [3] and [3] show typical angle- and energy-resolved 
photoelectron spectra for the s and the p x / 'p y -sta,tes. One 
clearly recognizes, for all cases, the formation of the ATI 
plateaus with distinct photon peaks separated by the 
photon energy u>. The isolated photon lines can be iden- 
tified for more than five orders of magnitude above the 
numerical noise level. In the case of the s-state two dif- 
ferent photon energies are plotted for illustration: Fig. [3] 
(b) shows the angular distribution for excitation with an 
energy of 27.2 eV and Fig. [3] (a) the distribution of pho- 
toelectrons ionized with 13.6 eV radiation. The energetic 
cut-off in both spectra is approximately the same, which 
can be attributed to the corresponding small pondero- 
motive potential in both cases. This corresponds for the 
case a) of smaller photon energy, lu = 0.5 a.u., to twice 
as many photon peaks, compared to b). 

A first inspection reveals that the photoelectron is 
ejected from the s-state predominantly along the polar- 
ization axis of the exciting field with a perfect angular 
symmetry around ip — and ip — n. In the case of 
both p-states one finds a small deviation from this sym- 
metry, best seen in Fig. [4] (a) and (b) at high photon 
orders whereas the main distribution — at least for high 
photon orders — is also ejected in directions near to the 
polarization axis of the electrical field. 

For further investigation of the processes involved, it is 
useful to construct the angular dependence of each sepa- 
rate photon line from the full spectrum obtained numer- 
ically. For this purpose we integrate the energy-resolved 
spectrum over an energy interval E . . . E+AE where AE 
is chosen according to the width of the photon line consid- 
ered. Since no channel switching or high-lying Rydberg 
resonances are accessed (due to the direct ionization with 
one single photon in the lowest photon order) the sub- 
structures in the photon peaks are not of importance in 
our case and are integrated out. This procedure is per- 
formed separately for each photon energy and intensity 
for every occuring photon peak in the spectrum. 

The angular dependence of photoelectrons initially 
bound in the isotropic 2s-state is shown in Fig. [5] for dif- 
ferent photon energies and photon orders. The columns 
correspond to the absorption of a fixed number of pho- 
tons Nph whereas the photon energy is varied for each 
row from lu — 2.5 a.u. (top row) to lu — 0.3 a.u. Each cho- 
sen energy is sufficient to ionize the electron directly by 
absorption of one single photon, cf. Fig. [2j One clearly 



sees the increased absorption of additional photons in 
the continuum with decrease of the photon energy, as 
discussed above. For the case of lu — 2.5 a.u. only four 
photon peaks can be identified in the spectrum whereas 
for lu = 0.3 a.u. more than n — 15 single photon peaks are 
found. The intensity of the laser pulse was chosen to be 
/ = 3.5 • 10 14 W/cm 2 for all photon energies. Therefore 
the ponderomotive potential U p = Eq /4lu 2 increases from 
top (U p = 0.0004 a.u.) to bottom (U p = 0.028 a.u.) but 
it is small compared to the binding energy and the pho- 
ton energy of the laser field. Therefore, all higher-order 
ionization processes (Nph > 1) are due to multiphoton 
absorption. Tunneling effects are expected to play only 
a negligible role. 



B. Angular dependence of different photon lines 
for p-states 

To investigate states with different angular momen- 
tum, as e.g., the 2p-states with I = 1, one has to keep 
in mind the orientation of the orbital with respect to 
the laser polarization axis. To construct the angular 
dependence of the photoelectrons for randomly aligned 
states it is necessary to address the specific orientational 
effects a single orbital shows in its PADs. In order 
to use the methods developed above we combine the 
photoelectron spectra of different states, i.e. for our 
special case (of the p-subspace) both 2p x and 2p y -states. 
Since the imaginary time stepping method constructs 
an orthogonal set of eigenstates, which are arbitrarily 
aligned with respect to the laser field, the dependence 
of the PAD on the explicit orientation of the 2p x and 
2p y states has to be investigated. For this purpose we 
consider, in addition, two rotated, orthogonal states 
2p' x and 2p' y (cf. Fig. [7j). Technically, the rotation 
of the states obtained numerically is performed by 
superposition of both original states, 2p x and 2p y , in 
such way that the final orientation of the new states is 
rotated by 45° compared to the original orientation. 

The obtained PADs for both cases are given in Fig. [7] 
The combined intensity of both states, I(p' x ) + I(p' y ), 
shows exactly the same behavior as the intensity of the 
non-rotated states I(p x ) + I(p y )- Therefore, special 
orientation-dependent features in the total angular 
distribution of photoelectrons can be excluded. 
Knowing this, it is sufficient to calculate the PADs of 
the p x and p y states. Since for atoms in the gas phase 
the orbitals occur randomly oriented, for every given 
state its orthogonal complement can also be found. We 
assume here, that the electrical field of the exciting 
laser is changing rapidly enough (high photon energies), 
such that the dipole moment of the atomic orbital 
does not play an important role and the atom is not 
aligned in some special way with respect to the laser 
polarization axis. Because the total PAD of each pair 
has the same shape, it is only necessary to know the 
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FIG. 5: Angular dependence of photoelectrons for different photon orders Nph ^columns) and different photon energies ui 
(rows, given in atomic units) for the 2s-state at an intensity of / = 3.5 • 10 14 W/cm . The laser field is polarized in ^-direction 
(horizontally) . 



ionization behavior of one single pair. In the following, 
we will therefore consider only two orthogonal states for 
the investigation of the total angular distributions of 
photoelectrons for p-states. 

Fig. [6] shows the obtained PADs from our TDSE sim- 
ulation in dependence on the photon energy (rows) and 
the photon order (column) for the 2p-state, similar to the 
case of the 2s-state, Fig. [5] But for this discussed 
above, each angular distribution contains three different 
graphs. The oriented p^-state (blue line), the oriented p y - 
state (red line) and the total intensity, plotted as an enve- 
lope (black line). As in the 2s case, one easily recognizes 



the emergence of additional photon lines with decreas- 
ing photon energy and the occurence of richly structured 
PADs. 

Additionally, the appearing of three regimes corre- 
sponding to different orientation of the PADs with re- 
spect to the laser polarization axis can be pointed out: 
in the first photon order (Nph = 1) the maximum inten- 
sity for ejected photoelectrons is oriented perpendicular 
to the laser polarization axis (left-most column in Fig.|6|. 
In the third regime (large numbers of absorbed photons) 
the PAD is aligned along the field polarization axis (right 
columns in Fig. [fjj. Finally, in the intermediate region, 
a transition-like behavior is observed which is character- 
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FIG. 6: (color) Angular dependence of different photon orders Nph (columns) for different photon energies u (rows, given in 
atomic units) investigated for the 2p-states at a laser intensity of I = 3.5 • 10 14 W/cm 2 . The colored lines indicate the angular 
distribution of photoelectrons of the single 2p x (red lines) and 2p y (blue lines) state. The laser field is polarized in a;-direction 
(horizontally) . 



ized by complex angular modulations. 



C. Intensity dependence of PADs 

Since scattering effects should become more prominent 
at higher intensities, or vice versa, decrease for low inten- 
sities, we calculated the PADs for different intensities of 



the laser field. The results are shown in Fig.[8]for a fixed 
photon energy lu — 1.69 a.u.. The columns again indi- 
cate the absorption of single photons and the rows mark 
different intensities from low (top, / = 3.5- 10 10 W/cm 2 ) 
to high (bottom, / = 3.5- 10 16 W/cm 2 ) [3T]. Again, with 
variation of the intensity, the number of observed pho- 
ton orders in the spectrum varies: from linear excitation, 
by absorption of one single photon with perfect dipolc- 
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FIG. 7: Comparison of the photoelectron angular distribu- 
tions for different orientations of the initial states. The upper 
figures give the densities of the rotated states (a) 2p' x and (b) 
2p' y . The PADs of the non-rotated 2p-states, cf. Fig. [2] are 
shown in (c) and for the rotated states in (d). The thinTines 
indicate the data for each individual orientation whereas the 
bold line displays the total intensity I(2p x ) + I(2p y ) in (c) 
and I(2p' x ) + I{2p' y ) in (d), respectively. 



like shape (top-most row in Fig. [8| in the perturbative 
regime, to the absorption of many photons with a richly 
structured angular dependence of the photoelectrons for 
different photon orders (lower rows in Fig, [8]), indicating 
scattering effects of the ejected electron with its parent 
ion. This trend is typical for all initial states (2s and 2p 
states). 



D. Physical explanation of the PADs in terms of 
scattering processes 

First of all, the shapes of the angular distributions have 
a qualitatively similar structure as found in experiments 
for IR photon energies [13] ■ But for all cases considered 
within the present work (except for very high-intensity 
excitations shown in Fig. [8j lowest row) the well-known 
ATI plateaus with characteristic energies of 2.5U p , 4.5U p 
and 10£/p are energetically located below the first ob- 
served photon peak, since the ponderomotive energy is 
small, U p < 0.03 a.u.. Therefore, the structures in the 
angular distribution cannot be attributed to the same ef- 
fects as in the IR photon case, i.e., to the onset of the 
ATI-plateau HHE2]. 

Nevertheless, scattering mechanisms should play 



the dominant role in the formation of the patterns 
observed. This we conclude from the fact that the 
richly structured shapes of the PADs, Figs. [5] [6] and 
[8] disappear both, for very large photon energies, and 
for simultaneous absorption of many photons at low 
intensities, respectively, where the electron is lifted high 
into the continuum and the absorbed field energy is 
converted directly into translational motion (kinetic 
energy) away from the atom. In this case, it is unlikely 
that the electric field, which changes very rapidly, drives 
the electron back to undergo a rescattering process. An 
analogous argument holds for the case of low intensities. 
The electrical field of the laser is not strong enough to 
modify the path of the fast-travelling electronic wave 
packet. In the following we will examine the obtained 
PADs more in detail to isolate traces of such scattering 
events. 

The shape of the PAD for the multiphoton ionization 
of order n is often fitted to a sum of even Legendre poly- 
nomials P21 [33] 



4?r 



1 + 2$? p 2i(cos9) 



(4) 



/=i 



to correctly describe the angular momentum the outgoing 
electron wave carries where 9 is the angle between the 
electron and the field polarization axis. In the case of 
the absorption of one single photon Eq. Q transforms 
into the well-known formula of Cooper-Zare [34] , 



(3P 2 (cos6)} 



(5) 



with the dipole-anisotropy parameter (3 and the angle- 
integrated ionization cross section a. Eq. (|5| describes 
the PAD for the ionization of randomly oriented atomic 
or molecular systems by a linearly polarized laser field. 
Its shape is, therefore, completely described by the two 
quantities (3 and a. 

If we now assume that, in the other regimes, scattering 
effects, and therefore changes of the angular momentum 
of the electron, influence the PADs, additional terms in 
Eq. Q beyond the first order, Eq. are needed. To 
this end, we fit a sum of Legendre polynomials P21 with 
I up to four via a least square fitting routine to our sim- 
ulation data. 

We will concentrate our discussion on the 2s state. In 
principle, the same results and arguments hold for the 
2p x - and 2f> J/ -states. The fitting results, i. e., the first 
four fit parameters 02.. .8, in dependence on the photon 
energy and the intensity for the first two photon lines 
are given in Fig. [9] and Fig. |10[ Higher-order processes 
would involve additional polynomials in Eq. |4) making 
the fitting procedure less reliable. 

First, we will discuss the behavior for the limit of large 
photon energies or low intensities. This is the regime 
where scattering effects are expected to play only a minor 
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FIG. 8: Photoelectron angular dependence for different photon orders Nph (columns) and different intensities (rows, electrical 
field amplitude Eq given in atomic units) for the 2s-state at a fixed photon energy of u> = 1.69 a.u. The laser field is polarized 
along the s-direction (horizontally). 
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FIG. 9: (color online) Parameters fii of Eq. ^ calculated via 
least square fit to our numerical data for the first (b) resp. 
second (a) photon peak of the PADs of the 2s-state vs. the 
photon energy u) for a fixed intensity of / = 3.5- 10 10 W/cm 2 . 



FIG. 10: (color online) Intensity dependence of j3i parameters 
in Eq. Q for the 2s-state at a fixed photon energy of lu = 
1.69 a.u. Again (b) gives the result for the first photon peak 
(single-photon absorption) and (a) the data for the second 
photon line. 



role in the angular distributions of photoelectrons. And 
indeed, the intensity distribution follows the Cooper-Zare 
formula, Eq. (j5j, in the case of single- photon absorption 
[Fig. [9](b) and Fig. 10 (b)] with an anisotropy parameter 
of = 2.0. This corresponds exactly to the theoretical 



value for dipolc radiation. All higher-order parameters 
0i>2 vanish. Analogously, the first two parameters 02 
and 04 are required to describe the distribution in the 
case of the second photon line, as it is expected from 
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N Ph = l 



N Ph = 2 



uj = 2.5 



u> = 1.0 



0.7 





order processes become relevant. This observation is 
condensed in the fact that higher-order polynomials in 
Eq. Q are accessed (/?4 and (3 e in Fig. 



10) 



FIG. 11: (color online) Comparison of PAD fitted according 
to Eq. (I| (black solid lines), with the accurate numerical data 
([red] dashed lines) for some selected typical cases of u>, see 
Fig. 



Eq. (|4| for n = 2. 

If the photon energy is decreased, the initial dipolc 
shape gets lost and the cocffccicnts of higher-order poly- 
nomials, 04, fie, . . . , increase where the next higher-order 
contribution is filled one after another. This behavior be- 
comes most obvious in the case of the first photon line, 
cf. Fig. [9] (b), where the /?2 parameter decreases from 
its dipole value #2 = 2.0, and the magnitude of the term 
/?4 is becoming larger (uj rs 2.5 ... 1.5 a.u.). In the fol- 
lowing step (uj < 1.5 a.u.), the third-order parameter (3§ 
gets involved to describe the PAD correctly. This behav- 
ior can be spotted directly in the shape of the PADs, cf. 
Fig. [5] first column Nph — 1, where up to a value of 
uj « 1.5 a.u. the dipole shape dominates and only small 
deviations are present. With decreasing photon energy 
the PADs become compressed along the x-direction and 
some parts of the electrons are scattered in different di- 
rections giving at low photon energies finally rise to the 
observed complex structured PADs. 

Accordingly, due to the higher kinetic energy of 
the photoelectrons being ionized by the simultaneous 
absorption of two photons and the corresponding smaller 
possibility of rescattering, the occupation of higher-order 
terms in Eq. Q sets in at smaller photon energies 
as can be seen in direct comparison of the (red) lines 
for 02 in Fig. [9] (b) and (a). For even lower photon 
energies, uj < 0.7 a.u. the angular distribution becomes 
very irregular, cf. Fig. [5] lowest rows, and the fitting 
procedure fails. 

The same behavior can be observed by increasing the 
intensity of the laser field, cf. Fig. [lOj where also the 
higher-order processes are accessed one after another. As 
in the previously mentioned case of variation of the pho- 
ton energy, the modifications can be directly spotted in 
the PADs, cf. Fig. [8] Again, the dipole-like shape at low 
intensities is, with increase of the laser intensity, com- 
pressed along the polarization axis and parts of higher- 



This lets us conclude, that scattering processes and, 
therefore, a modification of the angular momentum of 
the outgoing electron with the corresponding different 
angular dependencies, are dominant and responsible for 
the complex photoelectron angular distributions. The 
richly structured PADs are, therefore, also in the case 
of high photon energies at high intensities (but small 
ponderomotive forces), a consequence of the scattering 
processes the electron undergoes on the ion on its way 
from its creation by (multi) photon ionization to the 
detector. 

A final remark shall be made on the success of the fit- 
ting procedure, cf. Fig. [TT] Given the complex structure 
of the PADs, it is remarkable that the fit of only four Leg- 
endre polynomials is sufficient for resolving the main con- 
tributions to the PAD, at least for photon energies larger 
than ui = 0.5 a.u. For large photon energies (uj — 2.5 a.u.) 
the shape of the PAD is completely describable by such 
a series, the (red) dashed line of the numerical data and 
the black solid line of the fitted function show almost no 
deviations for both photon orders under investigation. 



In the transition- like region, in Fig. 11 represented by 
the case of uj = 1.0 a.u., first small substructures can 
be spotted, being not captured by the fitted polynomial. 
But the main contributions are still resolved within high 
accuracy. In the last case of uj — 0.7 a.u. given, the for- 
mation of jet-like structures at angles of approximately 
60° off the polarization axis of the laser plays against 
the fitting procedure, but still, the main contributions 
are accounted for. For even smaller photon energies, the 
strongly peaked, very sharp jet-like structures, cf. Fig. [5] 
require to extend the expansion, Eq. Q to orders higher 
than / = 4. 



V. CONCLUSIONS AND OUTLOOK 

In this paper we have investigated the angular distribu- 
tions of photoelectrons being excited by linearly polarized 
laser fields at high photon energies which are (or will soon 
be) available at free electron laser facilities. The observed 
PADs are, depending on the photon energy uj and the 
electrical field E of the laser, richly structured showing 
prominent side lobes and jets which were observed before 
in the IR regime. But, due to the small ponderomotive 
energy, the mechanisms cannot be attributed to the well- 
known formation of ATI plateaus. Nevertheless we have 
provided clear evidence that these effects are caused by 
scattering of the ionized electron on the ion when it is 
driven back by the laser field. This has been shown in 
terms of the analytical fit formula, Eq. Q , where subse- 
quently higher-order terms are needed as the possibility 
of scattering events increases. 

While our results where obtained for a 2D model atom, 
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we expect that the main features will survive in the 3D 
case. The present results are of relevance for single atoms 
and molecules (see e.g. [35]). They also apply to atom 
ensembles because, as we have shown, the total PAD of 
electrons ionized from all p-states is independent of the 
orientation of the atoms relative to the laser field. The 
high sensitivity of the PAD to the field strength allows 
for accurate detection of the laser intensity at the place of 
the atoms. Furthermore, the PAD in the UV and XUV 
regimes will be of high interest to pump-probe experi- 
ments with femtosecond to sub-femtosecond time reso- 
lution, e.g. [35] 13Z]. Finally, it will be of high interest 
to extend the present analysis to many-electron atoms 



by including correlation effects. This can be done by 
using e.g. quantum kinetic equations or non-equilibrium 
Green's function methods [35] [3Jjl 001 [H] which is subject 
of ongoing work. 
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